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ABSTRACT 



Although the internal linear combination method (ILC) is a technique widely used for the separation of the Cosmic 
Microwave Background signal from the Galactic foregrounds, its characteristics are not yet well defined. This can lead 
to misleading conclusions about the actual potentialities and limits of such approach in real applications. Here we 
discuss briefly some facts about ILC that to our knowledge are not fully worked out in literature and yet have deep 
effects in the interpretation of the results. 
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1. Introduction 

A widely used approach for the separation of the cosmic mi- 
crowave background (CMB) from the diffuse Galactic back- 
ground is the internal linear combination method (ILC). 
For instance, this method was adopted in the reduction of 
the data from the Wilkinson Microwave A nisotropy Probe 
(WM AP) satellite for CMB observations (|Bennett et al. I 
|2003| ). Its success is due to the fact that, among the sep- 
aration techniques, ILC calls for the smallest number of a 
priori assumptions. If the data are in the form of No maps, 
taken at different frequencies and containing Np pixels each, 
the model on which ILC is based is 



solution of type 



Na 



(2) 



If the constraint Wi — 1 is imposed, Eq. ^ becomes 



No 



(3) 



5«(p)-6c(p) + ^f (p)+A-«(p). 



(1) 



Here, S^^\p) provides the value of the pth pixel for a map 
obtained at channel " i " Q, whereas Sc(p), Sf^\p) and 
J\f^^^{p) are the contributions due to the CMB, the dif- 
fuse Galactic foreground and the experimental noise, re- 
spectively. Although not necessary, often it is assumed that 
all of these contributions are representable by means of sta- 
tionary random fields. Moreover, without loss of generality, 
for ease of notation the random fields are supposed as the 
realization of zero-mean spatial processes. The basic idea 
behind model H]) is that, contrary to the components that 
form the Galactic background, CMB is independent of the 
observing channel. ILC exploits this fact averaging No im- 
ages {S^'^^p)}'^^l and giving a specific weight Wi to each 
of them so as to minimize the impact of the foreground 
and noise ([Bennett et al. Il2003f) . This means to look for a 



Now, from this equation it is clear that, for a given pixel 
"p" , the only variable terms are in the sum. Hence, under 
the assumption of independence of &cip) from Sf^\p) and 
J\f'^^^{p), the weights {wi} have to minimize the variance of 
&cip), i.e. 

{wi} = argmin 



VAR [6c(p)] +VAR 



J2Msl'\p)+Af^'\p)) 



(4) 



where VAR[s(p)] is the expected variance of s{p). 

If S''*-' denotes a row vector such as S''*-' = 
[5'W(l)^5'W(2),...,S'W(iVp)] and the No x Np matrix S 
is defined as 

/ 5(1) \ 



In the present work, p indexes pixels in the classic spatial 
domain. However, the same formalism applies if other domains 
are considered, for example, the Fourier one. 



then Eq. ^ becomes 



S = &C + S{+J\f. 



(5) 



(6) 
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In this case, the weights are given by (|Eriksen et aL |[2004[ ) 

(7) 



where Cs is the No x No cross-covariance matrix of the 
random processes that generate S, i.e. 



Cs = nss^], 



(8) 



and 1 = (1, 1, ... , 1)"^ is a column vector of all ones. Here, 
E[.] denotes the expectation operator. Hence, the ILC esti- 
mator takes the form 

©c = w^S, (9) 
= al'^Cg^S, (10) 

with l^w = 1 and the scalar quantity a given by 

a=[l^C-^'l]-\ (11) 

In practical applications, matrix Cs is unknown and 
has to be estimated from the data. Typically, this is done 
by means of the estimator 



1 



-ss' 



(12) 



In this case, the ILC estimator is given by Eqs.([9|l- (fTT|) with 
Cs and w replaced, respectively, by Cs and 



Cs 1 
l^Cs 1 



(13) 



2. Some unfocussed points about ILC 

In spite of its popularity, various questions concerning ILC 
appear not yet well fixed. A first issue is linked to the 
fact that the estimation of the power-spectrum of CMB 
is a problem deeply different from the separation of this 
component from the Galactic foreground. Indeed, the esti- 
mation of the second-order properties of a stochastic sig- 
nal, though contaminated by noise, is an easier task than 
its recovery. Here, the point is that, if one is interested 
in the spatial distribution of the CMB emission, then in 

Eq. ([T]) the term S^^^ (p) + 7V^*^ (p) cannot be considered 
as a single noi se component. Often such assumption is 



made (e.g., see [E riksen et al. 112004: iHinshaw et a"]~ll2007t 
iDelabrouille et al . 2009: Dick et al. Il2009[ ) since in this way 



the problem is reduced to the separation of two components 
only and no a priori information on the "global" noise is 
required. Actually, this procedure can lead to wrong con- 
clusions. For example, since all of the component in the 
mixtures S are supposed to be the realization of a zero- 
mean random processes, from Eq. ^ one could derive that 

E[©c|©c] = ©c + w'^E[Sf] + w'^EiAf] = ©c, (14) 

i.e. the ILC^stimator is unbiased0. This is not correct: the 
claim that ©c is unbiased requires to prove that 



E[©e|©c, Si\ = ©c + W^Si -f W^nN\ = ©c 



(15) 



The reason is that Si is a fixed realization of a random 
process. There is no way to obtain another one. Even if 
observed many times (under the same experimental condi- 
tions) the foreground components (for instance the Galaxy) 
will always appear the same. Only the noise component M 
will change. Indeed, apart from very specific situations, in 
general I LC has to be expected to provide biased estimates 
of CMB (jVio fc Andreani 2008 1. Here, we stress that the 
biased results provided by ILC is a well established fact 
only relati vely to the estima tion of the power-spectrum of 
CMB Ce.g. lSaha et.~ani2008D . 

A second issue, often neglected in literature, is the be- 
lief that, in order to work, ILC does not require any hy- 
pothesis about S j . Aga in, this is not correct. As proved by 
I Vio fc Andreani I (|2008fl . even in the case of noise-free data, 
one may hope to obtain an unbiased estimate of the CMB 
only if Si is given by a linear mixture of the contribution 



of No physical processes {©j}^i, i.e 



s\ 



(16) 



with a.y constant coefficients and No < Nq. This means 
that for the jth physical process a template &j is assumed 
to exists which is independent of the specific channel " i " . 
Moreover, the number of these channels has to be greater 
than that of the physical processes. Inserting Eq. into 
Eq. ^ one obtains 



S = A&+Ar 



with 



© 



/©c 

©1 

©2 



(17) 



(18) 



and 



/I 



A = 



Vi 



an ai2 



0.21 0,22 



ONal aw 2 



02N^ 



ON„N, 



(19) 



This is the only model that ILC is able to handleEl- This fact 
should suggest extreme caution when dealing with ques- 
tions as the nonCaussianity of CMB. Indeed, the results 
produced by ILC are extremely sensitive to the presence of 
components that are not modelizable by means of Eqs. (|17p - 
(|19p . By instance, without renouncing the cosmic origin of 
C MB, this could exp lain the spatial associations claimed 
by IVerschuurl ()2007D between the interstellar neutral hy- 
drogen (HI) emission morphology and small-scale structure 
observed by the Wilkinson Microwave Anisotropy Probe 
(WMAP). 

A third issue is that, using the weights w, it is im- 
plicitly assumed that Cs is a well conditioned matrix. 
However, there are various practical situations where this 
condition is not satisfied. By instance, this can be expected 
for high signal-to-noise ratio observations at high Galactic 
latitude, i.e. when noise Af is negligible and CMB is by far 



2 Here, E[a\b] = [E[a(l)|6(l)], E[a(2)|6(2)], . . . , E[aiNp)\b{Np)]f 
with E[a(p)|b(p)] the conditional expectation of a{p) given b{p). 



^ If data are contaminated by noise, i.e. A/" 7^ 0, then in Eq. (O 
Cs has to be substituted with Cs — Cm ■ 
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the dominant component with <^ I, i — 1,2^ No, 
j — 1,2, ... , Nc. Indeed, since 

Cs = ACeA^, (20) 

when A is badly conditioned it has to be expected that the 
same holds for Cs- As a consequence, an imprecise cali- 
bration of the No maps, in such a way that the true value 
of the entries in the first column of A are different from 
"1" , will be amplified with cata strophic consequences. This 
is what lead iDick et al. I (|2009l ) to the conclusion that, in 
the case of high signal-to-noise ratio observations, ILC is ex- 
tremely sensitive to calibration errors. If true, this should be 
a quite troublesome situation since it limits the usefulness 
of the observations in a high signal-to-noise ratio regime. 
Actually, this problem can be easily avoided if in Eq. ([7]) 
Cg^ is substituted with the corresponding Moore-Penrose 

inverse Cg. These arguments hold also when the weights w 

— 1 

are used and Cg substitutes Cg in Eq. (|13p . This is clearly 
visible in Figs. [1] and [5] that show the results provided by 
ILC in the case of a simulated high Galactic latitude ob- 
servation at three different frequencies, say 30, 44, and 70 
GHz. The region analyzed is a square patch (400 x 400 pix- 
els) with side of about 24°, centered at / = 90°, 6 = 45° 
(Galactic coordinates). Three components have been con- 
sidered, i.e. CMB, synchrotro n and dust. More details can 
be found in IVio et al.1 (|2003[ ). The effects of the different 
point spread function for the various frequencies have been 
neglected. The top panels of Fig. [T] shows the templates rel- 
ative to the phsyical processes that have been used to form 
the three linear mixtures shown in the bottom panel of the 
same figure. No noise has been added. Perfect frequency 
scaling has been assumed, i.e. model (|17p - (fT9)) holds, but 
a calibration error of 10% has been imposed in all of the 
channels. The assumption of perfect frequency scaling con- 
stitutes only an approximation. However, since CMB by far 
constitutes the dominant component, this is of secondary 
importance. In fact, the condition number of the resulting 
Cs is 1.6 X 10^, i.e. this matrix is ill-conditioned. The re- 
sults obtainable by ILC using the original weights and 

their version computed using are shown in Figj^l It is 
evident that, if not properly addressed, the ill-conditioning 
of Cs has catastrophic consequences on the quality of the 
separation. 

3. Conclusions 

From the analysis presented above, it is evident that the 
results provided by ILC have to be interpreted with extreme 
caution. This techniques suffers many drawbacks that, if 
not properly taken into account, can lead to misleading if 
not wrong conclusions. In particular, ILC should be used 
only in situations where matrix Cs is (close to be) singular, 
i.e. only when one can be certain that model (|17p - (fT51) holds 
with No > Nc + 1. In the contrary case, the separation 
operated through ILC has to be expected inaccurate. 
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Fig. 1. Original templates (top panels) and the corresponding three mixtures (bottom panels) used in the numerical 
experiments described in the text. 
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Fig. 2. ILC solutions obtained using, respectively, the Moore-Penrose inverse of Cs (top-left panel) and the original 
weights (fT3|) (bottom- left panel). For reference, the target CMB component is also presented (top-right panel). 



